WHAT IS CLAIMED IS: 

1. A method of approximating an FIR filter with low-order 
linear-phase IIR filters by the rational Arnoldi algorithm with adaptive orders 
containing the following steps: a. initialize the first vector of the Krylov 

5 sequence for each expansion point; 

b. in the yth iteration of the algorithm, choosing an expansion 
frequency such that the frequency gives the greatest difference between the 
(/+l)st-order output moment of the original FIR filter H(z) and that of the 
lower-order IIR filter H(z) ;c. after the choosing the expansion point in j th 
10 iteration being determined, the single-point Arnoldi method applied at the 
expansion point to generate the new orthnormal vector; and 

d. determine a new residual at each expansion point for next 

iteration; 

whereby, after the giving total iteration number of the algorithm, 
15 outputting the resulting orthogonal projection matrix. 

2. The method as claimed in claim 1, wherein the exact expression of 
output moment errors between the ] t th-order moments H Ui \z t ) and 
H^{zi) at each expansion point z t are expressed as follows: 

| H (A ) (z, ) - H < A > ( Zi )HKc T r & -» (z f ) | , 

20 where 

K( z i) ^riylK^ 1 ^ 2 ^! * s t " ie norma li zat i° n coefficient when an expansion 
frequency z t is selected in the jth iteration; 

vector c contains the last n impulse response coefficients of a FIR filter with 

18 



length and 

r^~ l \zi) is the residual vector in the G"l) st iteration of the disclosed 
adaptive rational Arnoldi algorithm at the expansion frequency z t . 

3. The method as claimed in claim 1, wherein the heuristics of 
selecting expansion frequencies in advance for the proposed rational Arnoldi 
method are given by 

(a) low-pass filters: the proposed method with the expansion point 
^i=0; 

(b) high-pass filters: the special structures of state-space matrices used 
to present the duality between low-pass and high-pass filters; let 
state matrices become A=-A, b = b , c = c , and h 0 = -h 0 , the 
expansion point co x = 0 chosen to perform the Arnoldi algorithm; 
when the corresponding orthonormal matrix V q is obtained and 

then the high-pass IIR filter, which satisfies the same specifications 
as the original FIR filter; and 

(c) band-pass/band-stop filters: the passband edge and the stopband 
edge frequenciesbeing the appropriate candidate expansion points in 
meeting the specifications of the design, and other expansion points 
with uniform spacing being recommend to be selected. 
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